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Abstract 

We derive the tunneling rate for paramagnetic molecules in the context of a collective 
spin model. By means of path integral methods an analytical expression is derived. 
Given the very large spins in question (s ~ 3000^), the observation of magnetization 
changes due to pure unitary tunnel effects is unlikely. 



1 Introduction 



The spin tunneling in magnetic materials has atracted much attention since 
experimental evidences of this phenomenon were observed. Although tunnel- 
ing is a purely quantum effect, it is interesting to note how it still mani- 
fests in superparamagnetic particles of macroscopic size (of the order of tenths 
of nanometers). These "particles" are generally magnetic grains or clusters, 
whose spins are highly correlated, and they behave like they were a single 
large spin particle [1,2]. The Macroscopic Quantum Tunneling (MQT) is the 
tunneling of a single macroscopic variable (in the present case, the magneti- 
zation) through a barrier, between two minima of the effective potential of 
a macroscopic system. For small magnetic clusters, these minima correspond 
to two states of opposite magnetization (easy direction), and the barrier is 
proportional to the exchange anisotropy. When there is a repeated coherent 
tunneling betwenn the two wells, then we have a case of Macroscopic Quantum 
Coherence (MQC), and all spins have their relative orientation unchanged[3]. 
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The existence of tunneling removes the degenerescence of the energy levels 
of the system, and the tunneling rate is given by the energy difference be- 
tween these levels. A useful and well known technique for performing these 
calculations is the instanton[4\ , which is based on the Path Integral formalism 
developed by Feynman[5]. In fact, the application of the method, together 
with the Villain transformation to the tunneling of spins has been introduced 
by Enz and Schilling in ref [6] (see also the review by Schilling in [7]). This is 
a non-pertubative method that provides the energy of the fundamental states 
with great precision. The proposal of this work is to show how to apply the 
instanton technique to the problem of the spin tunneling of a superparamag- 
netic particle with a detailed discussion of all the necessary steps to find the 
energy sphtting of the lowest levels. 



2 The Hamiltonian of the system in terms of canonical operators 

We start here with a general Hamiltonian of a single large spin system with 
easy direction on the z-axis which has been used in the context of superpara- 
magnetic particles [8], given by 



with the anisotropy constants A and B satisfying A ^ B > 0. Sz and S± — 
Sx ± iSy are the standard spin operators. In the basis of the eigenvectors of 
and Sz, these spin operators satisfy 



To solve this problem, we shall use the Villain transform[9], which consists of 
defining the operators S+ and S_ in terms of and a new operator (f, 



H^-ASl + B{Sl + S^_), 



(1) 



Sz\s, m) — m|s, m) 



(2a) 
(2b) 



S±\s, m) = J{s ^ m){s ±m + l)\s,m± 1) 




(3a) 
(3b) 



with 



Sz] 



(4) 



The eigenvectors of {|0)}, are such that 



e' 



(5) 
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To expand 5*+ and S- wc need that ||5'^|| ^ s. This would be true, were the 
system's easy direction not on the 2;-axis, but on the y or x-axis. This direction 
is clearly arbitrary, so we choose it on the y-axis, and rewrite the Hamiltonian 
as 

H^-AS'^ + 2B{S',-Sl), (6) 

or equivalently 

H = j{Sl + Sl- S+S^ - S^S+) + 2B{Sl -^{Sl + Sl + S+S. + S.S+)). 

(7) 

Besides the commutation relation (4), we need another one, which can be 
derived from it and proves to be much useful, that is 

[S„e^'""^]^±me^''^'^ . (8) 

One implication of the above equation is 

/(5,)e±^^ = e±^^/(5,±l), (9) 

where f{Sz) is any function of S^. With these relations, we can write the 
Hamiltonian, in the "normal order" , in terms of these new operators, that is 



H = A[{s{s + 1) - 5^^)(cos2 ^ - 1) + iS, sin 2ip] + 
2B[Sl - {s{s + 1) - Si) cos^ if - iS, sin 20]. 



3 Transition amplitude 

We have to compute the following matrix element 

(0/,i/|0,,i,) = (0/|e-^(*/-*^)|0,). (11) 

This is achieved in a standard fashion by discretizing the time in intervals 
e = {tf — ti)/N and inserting different projection operators like/ d4>k\4'k){4'k\-i 
between each time interval. We get 

{(t>f^Mi^i) = Yira^j d0,j (0;|e-^|0^_i)...(0i|e-^|0,) (12) 
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Each one of the inner products above is given by 



ieH 



s c,irrik{'t>k-<t>k-l) 



(13) 



m,},=-s 



27r 



with 

H{mk,(f)k-i) 



\H\<Pk-i) 



{(pk\<Pk-i) 

= + 1) - ml) (cos^ - 1) + iruk sin 20fe_i] + 

2B[m1 — {s{s + 1) — ml) cos^ (pk-i ~ i'^^k sin 20fe-i] 

Grouping the mfe-terms and making the substitution 

{A-2B) sin20fc_i 



(14) 



TTT-fe + I 



2[A + 2B-{A- 2B) cos^ 



(15) 



we get, after transforming the summation into an integral and taking the hmit 
s ^ oo (see Eq. (2.13) of Johnson[10]), 



fe|e-^|0fc_i) 



27r 



gi(pfc+A)(0fc-</>fc_i) ^-ieH{pk,(t>k-i) 



(16) 



where 



H{pk, <f)k-i) 
M{4>k-i) = 



pI 



2M{cj>k-i) 



+ V{cj>,_,) 
1 



2[A + 2B-{A- 2B) cos^ 
V{(Pk-i) = {A- 2B)s{s + 1) cos^ - As{s + 1) 

and A = if s is integer or A = | if s is half-integer. Now we use Eq. (16) 
back into Eq. (12) and take the limit ^ oo to find 



{<i>f,tf\<f>i,u)^ lirn^j [n (n 



27r 



N 



X exp<^ie^(pfe + A) 



/ (pk - <Pk- 



= j V[(i)\D\p]ex^U j dt 



(p + A)0 



2M{(t)) 



H{pk, (pk-i 



V{cj>) 



(17) 
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Completing square for p^, and computing the gaussian integral, we get 



{4>f,tf\4>i,U) = j V[4>{t)]ewU J dt 



(18) 



Making a Wick rotation, we then have 



{<Pf\ e-^(-/--^) 10,) = J l)[0(r)] exp |- / dr 
= / mr)] 



3— 'S'b 



(19) 



from where we find the classical Hamiltonian and Lagrangian of the system 

H{p, (p) = [(25 - A) cos^ <f) + 2B + A]p^ + {A- 2B)s{s + 1) cos^ - As{s + 1) 



L(0,(/)) = 



0^ 



4[(2S - A) cos2 + 25 + ^1] 



(20a) 

(^ - 2B)s{s + 1) cos^ (f) + As{s + 1) 

(20b) 



4 Expansion of the action with variable mass 



Let us take s integer, so that A = 0, as it is appropriate for superparamagnetic 
particles. The case with A = | leads to a topological selection rule, already 
discussed in the literature[ll]. The Euclidean action is then 



Se = dr iM(0)02 + V{ct>). 



(21) 



The instanton is the solution which minimizes the Euclidean action. In the 
semi-classical limit, we make an expansion of the action around the instanton 
0, taking 



= + CiSifi 



(22) 



(repeated indices indicate summation). 5(pi are eigenfunctions of an operator 
X>op, and form an orthonormal basis with the boundary condition Sipi{±T/2) — 
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0. Expanding each one of the terms above, we obtain 



M{(f)) = M{4> + aScpi) = M(0) 



dM 



CiS(pi + 



2 (902 



1.2 



= (0 + CiS(pi) = + 2(t)Ci5(pi + CiSipiCj60j 



dV 

v{<P) = v{<p) + — 



CiOipi + - 



2 



Ci5(piCj6(pj 



CiS(piCj5(fj (23a) 

(23b) 
(23c) 



Substituting the expressions above into Eq. (21) and integrating by parts, the 
action becomes 



1 /•'^/2 

So + - dr ^ CiCj5(piVop5(pj 



The transition amphtude is then 

(0^|e-^(-/--^ 10,) = e-^"[det(I^op)]-'/'. 



(24) 



(25) 



A chain of instantons and anti-instantons is also a solution which minimizes 
the action. When we take these solutions into account (see [4]), we get finally 



n odd 



where T = Tf — ti and 



K — lim 



dct(-D, 



-1/2 



op; 



T-co Vdet(-a2 + a;2^ 



(26) 



(27) 



5 The instanton 



To find the instanton, we write the Euclidean action as 



Se^ j dT[-ip0 + //(p, 0) 



(28) 



Using simple variational calculus, we have that the 5p coeficient must vanish, 
and then 



i0 = 2[{2B - A) cos^ + 2S + A\p. 



(29) 
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The condition that the instanton is the solution of minimal energy gives 

H{pJ) = -As{s + 1) 

^ {2B - A) cos"^ 4) + 2B + A' ^ ' 

and so 



i4) = ±2V2B - A^s{s + 1) I cos 01 V (25 - A) cos^ (j) + 2B + A (31) 
We make k = |f^, such that — 1 < A; < 1, and then 

= ±2VA^ -AB^^s(s + l)\ cos 0|Vl + A;cos2 (32) 
Now, integrating the equation above, we have 

/ , "^'^ J = ± / 2^s(s + 1)VA2 -452 dr. (33) 

cos 1 + A; cos^ 

Making the substitution x — sin 0, dx = cos (p dcp, we can readily solve the 
integral to find 

dx - 1 In f + ^ ~ + A C34N 

(1 - x2)^l + k{l - x2) ~ 2 " V VI + A; - - x) 



It is easy to check that uj = yV"{(/)o) /M(0o) is equal to 2y^s(s + 1)VA^ - 4i?2 
(00 is the initial value of 0(r) equal to 0(±oo) = ±7r/2), and therefore the 
solution is given by 

~ . / VTTksmh[±uj{T-To)] \ 

(p — asm < ' ' = > (35) 

[yjl + {l + k) sinh'[±u;(T-To)]J 
^("^ ^^VlTfe cosher -ro)] ^''^ 



There are thus two instanton-like solutions, with positive (/+) and negative 
(/~) signal of Eq. (36), and the corresponding anti-instantons (/+ and /"), 
given by the positive or negative signal of the argument of the function sinh in 
Eq. (35). Figure 1 shows the difi^erent types of instantons and anti-instantons. 
They are necessary to the calculation of the Berry's phase later (section 6). 

Making the identification p — i^Js{s + 1) cos(6'), we obtain the classical angle 
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e 



, 2J-k/{l + k) , 
0{t) = acos ^^-4^ 1 . (37) 



Figure 2 shows the positive graphs of these functions. We can see clearly the 
expected result. For example, for the instanton represented in this picture 
(/+), we have (p = —7t/2 and 9 = 7t/2 when r — >• —oo, and (p = 7i/2 and 
9 — 7r/2 when r ^ oo, which are the states of minimal energy. The classical 
action can now be readily computed 



rT/2 i-k/2 . 

^0 — —i / dr(f)p — —i d(f)p 

J-T/2 J-tt/2 

So = ^Js{s + l)ln (^^ 



rT/2 i-k/2 
-T/2 J--n-/2 

^^^^ (38) 



5.1 Eingenvalue zero 

If we apply the operator X>op to the function 

we see the result is zero because the term in parantheses is simply the equation 
of motion of the system. Then we may conclude that is the eigenfunction of 
Pop with eigenvalue zero, and we must treat this eigenvalue separately. 

If (j) centered at tq is an instanton solution, then centered at tq + a must 
also be. Thus 

(pTo+air) = 4>to{t - a) ~ (pToir) - acp = 0ro(^) + cqSlpq, (40) 
with the normalization condition imposed by 

rT/2 j-T/2 -Jl 

/ dTbkpl = N \ dT^ = 1. (41) 
J-T/2 J-T/2 

The integration over Cq was already performed while computing 5*0 and there- 
fore must be taken off the determinant. However, we must correct the fact that 
the integration was over tq, and consequently we need a correct normalization 
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factor. A variation dro in the center of the instanton causes a variation 



A[0(r)] = dro. 

Since 0(r) = 0(r) + CiSipi, a variation in cq associated to 5ipQ, gives 

A[0(r)] = (ico^v^o- 
As 5lpq = l/v^0, then dco = dro, and 



K 




hm 



det'(r', 



-1/2 



op 7 



27r T^oo \^det(-92 + 0^2 



(42) 



The prime indicates that the eigenvalue zero must not be computed in the de- 
terminant and the factor l/^/27r comes from the definition of ©[^j = Yin dCn/ v^27r. 



5.2 Calculation of K 



From the theory of differential equations, the following relation[^ can be de- 
rived for the eigenvalue A = 0, 

det(P^p) ^ ij\T/2) 

where T>\^ are second order differential operators in the interval [— T/2,T/2] 
and ip^'"^ are eigenf unctions of these operatores, such that 

^^'2(-T/2) = 6 dr^^'\~T/2) = 1. (44) 

We multiply and divide Eq. (43) by Ao(T), which is the smallest eigenvalue 
of Pop and vanishes in the limit T oo. 4''^{r) = ^sinh[ci;(r + T/2)] is an 
eignefunction of the operator —d'^ + u"^ which obeys the conditions above and 
has the form lim^-^oo = e^'^ /2uj in the limit T — oo. So, in this limit, 

Eq. (43) becomes 

det(Pop) 1 ^ 2^^(T/2) 
Aei{-dl + u^)\o{T) Ao(T)e-^- ^ ' 



5.3 Calculation of iP{t) 



As already mentioned, is an eigenfunction of the operator Pop- Now we 
calculate this function explicitly starting from the instanton solution (Eq. 

2 A prove of this relation may be found in Coleman [4], Appendix 1 
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(35)). Taking the derivative with respect to r, we see that 

~, , uj\/l + k cosh (cur) 
"^^^^ " 1 + (1 + A;) sinh^ur)' 

In the hmit r — > ±00, which is the hmit of interest, this equation is simply 

lim 4)(t) = — ^^e-'^l^l 



and after normalization, we have 

7 . . 1 2a; 



NVTTk 



e-'^l^l = Ce-'^l"! = xi. (46) 



This function however doesn't fulfill the conditions imposed by Eq. (44), be- 
cause dr(pN{^ = —00) = 0. We then seek a solution yi, linearly independent 
of Xi, such that t/j — Myi + Nxi satisfy the boundary conditions. We write 
T^op — d'^+I{T)d+J'{T). If Xi and yi are eigenfunctions of X>op with eigenvalue 
zero, so 



y,[d^x, + I{T)dx, + J'{T)x,] =0 
x,[d^y, + I{T)dy, + J'{T)y^] =0 



Subtracting both equations 

d(yidxi - xidyi) + I(r)[yidxi - Xidyi] = 

and definig W{xi,yi) = yidxi — Xidyi, the Wronskian of the system and 
dF/F = I{t), we sec that 

dF 

dW{x,,y,) + —W{x,,y,)^Q ^^^^ 
d{FW) = ^ FW = constant. 

As "Dop = -Mdl - {drM)dr + J{t), we can identify F with M(0), and then 

M{(f))W{xi,yi) = constant. 

We normalize W{xi,yi) such that M{4>)W{xi, yi) = 2MquC'^, where Mq is the 
limiting value of M(0) when r — > ±(X). Thus, in this limit, W{xi, yi) = 2a;C^. 
Using the Wronskian, we get, in the limit r — > ±00, 

yi = ±Ce±-^ 

or equivalently 

yi = ^Ce'^l-l. (49) 
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We need now to find the coeficients M and N, in order that ijj satisfy the 
boundary conditions. They are 

g-wT/2 qU}T/2 

From the above, we finally have 



g-<^lr| 

|t| 



(50) 



and therefore we find the value of ip{T/2) = 1/ui. 
5.4 Calculation of Xq 

We start with the eigenvalue equation of the operator V^p, 
We write ipxir) in an integral form like 



Mr) = M) - D / ds[y,{T)x,{s) - x,{T)y^{s)]Ms) (51) 

J-T/2 

where D is such that the solution above satisfy the eigenvalue equation. Taking 
the first and second order derivatives of iP\{t), we have 

Mr) = V'(t) -D r ds[y[{r)x,{s) - x[{r)y,{s)]Ms) 

J-T/2 

- D [|/i(t)xi(t) - xi(r)yi(T)] ipxir) 
V ' 

=0 

Mr) = Mr) -D r ds[y';{r)x,{s) - x'l{r)yr{s)]Ms) 

J-T/2 

-D[y[{T)x^{T)-x[{T)y,{T)]Mr)■ 

^ V ' 

=W{xi,yi) 

Now, we substite it back into the eigenvalue equation, and use the equation 
of motion, -M^'{t) - {drM)^'{T) + J(r)^(r) = 0, to get 

[A - MDW{x^,y^)]Mr) = D ds!^M[y';{T)x^{s) - x';{T)y,{s)] 

+ {drM)[y[{r)x,{s) - x[{r)y^{s)] - J{r)[y,{r)x,{s) - x^{T)y,{s)]'^Ms) 

(52) 
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The term on right is equivalent to 

D f_^^^ ds\[-Mdl - {drM)dr + J{T)][yi{T)xi{s) - (r)7/i (s)] j^^A (s) (53) 

As Xi{t) and |/i(t) are eigenfunctions of V^p with eigenvalue zero, then the 
right side is null. Thus 

[X-MDW{x^,y^)]Mr)^0 
D = - - (54) 

Iterating Eq. (51), we have 

^;,(r) = iPir) -dJ ds[- ■ ■ ]ij{s) + J ds[- ■ ■] J ds'[- ■ ■ ]tP{s') + ■■■ (55) 
Since A ^ 0, then D, and thus 

Using the results for xi and yi, we have, for r = T/2 
1 A /-^/s 

= ^ - iJwc^ -^^I^"' ^?(^) - ' i-'*^'!- (^*" 

This integral can now be readly solved and gives a result equal to C'^/uje'^'^ . 
Therefore 

1 XpC/uje^T _ 
Ao = 4Moa;^e-'"^ (57) 

5.5 Value of K 



With the results obtained here so far, the value of K may be readily computed 

V 27r T-.00 \^det(-92 + u;^) y y 27r r-.oo \^ Ao(T) e'^^ J 

2^ 



Mo is given by lim<^^±^/2 (0) = [2(^ + 25)] ^. N is given by the nor- 
malization condition of Xi, and may be solved to give N — 4:Lj/{1 + k). 
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1 + A; = 45/(25 + A) and u = 2^s{s + 1)VA'^ - AB^. Hence the final re- 
sult for K is 



K=^=. (58) 



6 Berry's phase 



As seen in section 5, the model has two instanton-hke solutions and their two 
corresponding anti-instanton. Then, there are two possible paths to go from 
one minimum state to the other and two others to return. The action along 
any of the paths is the same up to a geometrical phase (Berry's phase) which 
is now computed. This phase is given by 

— exp is^l — exp d0 is| 



For the instanton /+ and the anti-instanton / , the Berry phase is 

f7r/2 
-7r/2 

and to the instanton I" and the anti-instanton 7+, 

-tt/2 



$5(/+) = $s(/") = exp I J^^^^ d0is| = e^'^" (59a) 



$5(7-) = $b(7+) = exp I d(/)is I = e"^'^" (59b) 



We must introduce these phases in the calculation of Eq. (26), because the 
contribution of each instanton is now multiplied by this phase. Taking account 
of the four possibilities of a 2-link chain (back and forth), the summation in 
Eq. (26) becomes 

(Ke-^° $B(7+)r)(7re-^o $s(7+)T) (T^e"^" ^b{I')T){K e'^" $b(7-)T) 

2 I ^ 2 

^ {Ke-^° $B(7+)r)(7re-^o ^b{I~)T) ^ {Ke-^° ^b{I')T){K e-^° $b(7+)T) 



2 2 

{Ke-^o T)2($5(7+) + <1.^(7-))($s(/+) + $b(/-)) 



(60) 



With the results obtianed in (59a) and (59b), we see that 

($b(7+) + $b(7-))($b(7+) + $b(7-)) = 4cos2(7rs) (61) 
and therefore in Eq. (26), we must replace e~'^°/^ by 2cos(7rs) e~^°/^ 
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7 Splitting of the ground state 



Now we can finally calculate the energy splitting of the ground state due to 
the ressonant tunelling between them. Summing up Eq. (26), we have 

from what we get the energies of the two lowest states to be 

E± = -a;±2/s:e^Ocos7rs 

Then the splitting in energy is 

AEi^t ^E+-E_^ 4/s:e-^o costts 

Replacing the values of uj and k in terms of A, B and s, and using that s ^ 1, 
such that s{s + 1) s + |, the final result will be 

(64) 



If we take the limit B/A<^1 and s » 1, the equation above is equal to 

Af^ fB_ 



A-Einst = ^ -r COSTTS. (65) 



This result may be compared with the expression obtained by Hartmann- 
Boutron[12] using other method and with the rigorous universal upper bound 
for the tunneling rate of large quamtum spins recently obtained [13] and inde- 
pendent of the form of the anisotropy, verifying an excellent agreement with 
both results. 

We can still test the result, comparing it with the exact value obtained via 
diagonalization of the Hamiltonian. In Table 1 we contrast these results, and 
it can be seen that, although the analytical calculation is valid for s ^ 1, even 
for small s, the discrepance between the values is small. 
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8 Conclusion 



We have found an analytical expression for the tunneling rate of a superpara- 
magnetic particle using the instanton method and compared it with others 
results obtained by different methods, and also with numerical calculations. 
The result proves to be in good agreement with all of them. It is worth noticing 
that for a typical superparamagnetic particle, with s ~ 3000 and / ~ 10~^, 
the tunneling rate is practically zero, and then, no tunneling should be ob- 
served in practice. However, the change in magnetization direction of such 
particles is observed, what suggests that we must take the interaction with 
the environment into account. Work along these lines is in progress. 
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Fig. 1. All possible instanton and anti-instanton solutions 




Table 1 

Comparison between the numerical and calculated results for the ground state en- 
ergy splitting for different values of s. Here we took / = 0,2. 



5 5,61125 xlO"^ 5,02226 x 10"'^ 

8 1,08198 x 10"^ 1,06013 x 10"^ 

10 1,48684 X 10"^ 1,51767 x 10"® 

12 1,93224 X 10"^° 2,04068 x 10"^° 
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